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The effect of oscillatory shear flows on turbulent transport of passive scalar fields is studied by 
numerical computations based on the results provided by E. Kim [Physics of Plasmas, 13, 022308, 
2006]. Turbulent diffusion is found to depend crucially on the competition between suppression due 
to shearing and enhancement due to resonances, depending on the characteristic time and length 
scales of shear flow and turbulence. Enhancements in transport occur for turbulence with finite 
memory time either due to Doppler and parametric resonances. Scalings of turbulence amplitude 
and transport are provided in different parameter spaces. The results suggest that oscillatory shear 
flows are not only less efficient in regulating turbulence, but also can enhance the value of turbulent 
diffusion, accelerating turbulent transport. 
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I. INTRODUCTION 

Turbulent transport is an important process, by which 
various physical quantities are rapidly mixed by the ad- 
vection of turbulent fluid. Examples include the mixing 
of chemical species, dissipation of magnetic fields, trans- 
port of angular momentum in astrophysical/laboratory 
plasmas, the mixing of pollution in the terrestrial atmo- 
sphere or even the transport of warm water in oceans. In 
fact, the effect of turbulent transport in our lives is fun- 
damental. Turbulent transport can however have an un- 
welcoming consequence, such as anomalous heat/energy 
loss in laboratory plasmas, leading to the degradation of 
plasma confinement. It is therefore crucial to understand 
the physics of turbulent transport, especially some means 
of controlling it. 

One of the promising mechanisms for quenching turbu- 
lence mixing is flow shear suppression P, 0, . The basic 
idea is that shearing of turbulent eddies by flow shear 
leads to a cascade of eddies into smaller and smaller scales 
down to the dissipative scale where the molecular dissi- 
pation efficiently removes the energy. That is, the shear 
flow rapidly generates small scales, thereby enhancing the 
overall dissipation of turbulent eddies which are responsi- 
ble for turbulent transport. As a result, both turbulence 
level and transport are reduced 0, 0, la, IS ■ I n f ac t> 
there has been accumulating experimental evidences for 
turbulence regulation by flow shear in laboratory plas- 
mas, which is now thought to be indispensable for the 
formation of transport barriers [H, 0] . A crucial question 
in understanding the formation of transport barrier in a 
variety of systems is thus, how much of turbulent trans- 
port is reduced by a given flow shear compared to the 
case in the absence of the flow shear. Quantitative theo- 
retical predictions for turbulent transport were provided 
in different types of turbulence models for steady 0, 0] 
or random shear flows 0, ■ 

In comparison with a steady shear flow, turbulence reg- 
ulation by time- varying shear flow is more complex, with 
its efficiency depending on characteristic time scales of 
the shear flow and turbulence such as their frequencies, 



decorrelation times, shearing rate, etc. 0, [S [T(|. For in- 
stance, if a shear flow oscillates rapidly with its frequency 
much larger than the turbulent decorrelation rate, its 
shearing is too incoherent to have any effect on turbu- 
lent transport. In contrast, if the flow shear is sufficiently 
strong such that its shearing rate involves the fastest time 
scale, it can be considered to be steady and thus has a 
similar effect on turbulent transport as a steady shear 
flow. For instance, in this limit, turbulent mixing of pas- 
sive scalar fields is reduced inversely proportional to the 
rms shearing rate Sl m as oc l/f2 m |9j. What happens 
in between these two extreme limits is far less obvious. 
While theoretical predictions in a scalar field model are 
provided by Kim [9(], they are given in terms of multi- 
ple integrals which could be reduced to simple analytical 
forms only in the extreme limits, yielding transparent 
scalings with fi TO only in these limits. Further quan- 
titative study is thus necessary in order to understand 
how efficiently an oscillatory shear flow controls turbu- 
lent mixing in general. 

It is important to emphasize that the aforementioned 
reduction in turbulent transport by flow shear results 
from the modification of the turbulence properties (e.g. 
enhanced dissipation of turbulence/fluctuations) when 
the flow shear can have no influence on the mean field 
profile, for instance, when a mean field, being uniform 
along the shear flow, varies only transverse to the shear 
flow. This is relevant for the study of turbulent transport 
where main the interest lies in the transport mediated by 
turbulence. In the case when a mean field varies along the 
shear flow as well as transverse to it, the shear flow has 
a direct influence on mean field since it can distort the 
distribution of the mean field, thereby rapidly generat- 
ing small scales which then get efficiently damped by the 
molecular diffusion [for example, see Refs. [U, HI, Ull]. 
That is, the same shear distortion that enhances the dis- 
sipation of turbulent eddies, leading to the reduction in 
turbulent transport, can speed up the diffusion of mean 
field if it directly operates on the mean field. Therefore, 
in general, when shear flow can modify the mean field as 
well as turbulence, the overall effect of the shear flow on 
turbulent transport will be determined by the competi- 
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tion between these two conflicting effects. 

It is also important to note that propagating oscilla- 
tory shear flows (waves) have often been invoked as a 
mechanism for transport [bj ]. Examples include inter- 
nal gravity waves which can transport momentum, mass, 
heat, etc in astrophysical and geophysical systems [e.g. 
see, Ref [IB]]. The transport by these waves is however 
a slow process compared to turbulent transport as the 
former requires non-ideal effect such as molecular dissi- 
pation. On the other hand, oscillatory shear flows can 
potentially mediate fast transport by destabilizing the 
equilibrium via parametric resonance [l6l . [l7| . In the 
presence of turbulence, the interaction between oscilla- 
tory shear flows and turbulent flows could lead to similar 
parametric resonance, as indicated by Kim Q, signifi- 
cantly contributing to the transport. Without this res- 
onance, the effect of oscillatory shear flows on laminar 
diffusion can be opposite since the shear enhanced mix- 
ing becomes negligible as the oscillation frequency of the 
shear flow increases above the rms shearing rate [l8j . 

The purpose of this paper is to perform a detailed 
quantitative study of the effect of oscillatory shear flow 
on turbulent transport of passive scalar fields. In par- 
ticular, we derive the asymptotic scalings of turbulence 
amplitude and transport with the rms shear strength via 
numerical integrations of the theoretical results derived 
by Kim [9( . We distinguish three different scaling regimes 
depending on the characteristic time scales of turbulence 
and oscillatory shear flow, identify the two types of res- 
onances (Doppler and parametric resonances), and then 
obtain asymptotic scalings with rms shearing rate valid 
in each regime by numerical computations. We show that 
turbulent diffusion depends crucially on the competition 
between suppression due to shearing and enhancement 
due to the resonances, thereby suggesting that oscillatory 
shear flows can enhance the value of turbulent diffusion, 
accelerating turbulent mixing. We note that our numer- 
ical integration has a great advantage of permitting a 
thorough parameter scan in parameter space, which is 
not easily accessible by direct numerical computations 
(e.g. in the limit of small dissipation). The remainder 
of the paper is organized as follows. In §11, we present 
the time-averaged turbulence amplitude and transport in 
dimensionless form. In §111, we identify the three scaling 
regimes depending on characteristic time scales by simple 
analytical examination. §IV contains the results of nu- 
merical computations and scalings of turbulence ampli- 
tude and transport with rms shearing rate. Discussions 
and conclusion are provided in §V. 



II. GOVERNING EQUATIONS 

We consider a passive scalar field model where a pas- 
sive scalar n is advected by a given turbulent flow u and 
shear flow Uo while being diffused by molecular dissipa- 
tion D. By quasi-linear analysis, the fluctuating scalar 
field n' evolves according to the following equation: 



{d t + U • V}n' = -v x d x N + DV 2 n'. 



(1) 



Here, iVo = (n) = Nq(x) is the large-scale component. 
We assume that the shear flow is in the y direction, vary- 
ing linearly in x, and is oscillatory in time of the form 
Uo = — xQ m sin (oj z t) y; tu z and £l m are the frequency 
and rms shearing rate of the oscillating shear flow. Note 
that since both Nq and Uo depend only on x while Uo 
is in the y direction, there is no direct effect of the shear 
flow on the mean field Nq (i.e. Uo ■ ViVo — 0). That is, 
the shear flow influences the mean field only indirectly 
through its effect on turbulence. 

The advection by a linear shear flow U (x, t) — —xil(t) 
results in the distortion of an eddy (i.e., wind-up), and its 
effect can be non-perturbatively captured by employing 
a time-dependent wavenumber k x {t) with the following 
transformation for n': 



n'(x, t) — n(k, t) exp {i(k x (t)x + k y y)} , 



(2) 



and similarly for v x , with k x satisfying an eikonal equa- 
tion 



8 t k x {t) = k y Vt{t) 



(3) 



It is worth noting that Eq. ^ clearly shows that the 
shear flow has no influence on the mean field (k y = 
mode) of our interest. The solution to Eq. (|TJ) can then 
be expressed as: 



h(k, t) 



-d x n 



f dt^hgik, t; ki, h)e- DQ ^h x {^ t($) 

J — oo 



Here, Q{t,ti) = f dt'[k^.(t') + k 2 ], and g is the Green's 
function given by 

g(k,t;ki,ti) = 6(ky-ki y )6 k x - k lx - k ly dt'Q(t') 

(5) 

The overall effect of enhanced dissipation due to shearing 
is embedded in the time integral of k 2 in DQ in Eq. 
(|3|). The detailed form of Q for Cl(t) = — tt m smuj z t was 
provided by Kim 0, to which the readers are referred. 
We here emphasize that Q grows at most linearly in time 
for the oscillatory shear flow as k x oscillates in time (see 
Eq. (O), with its shearing becoming effective only when 
it operates coherently before the oscillatory zonal flows 
change shearing direction, i.e. whenf2 m /w z > 1 [9(. This 
is distinctively different from the behaviour of Q in the 
case of mean shear flows or random zonal flows: (i) in 
the case of mean shear flows [J], k x grows linearly in 
time with Q cx t 3 ; (ii) in the case of random zonal flows 
0] k x can be considered as a random process with the 
time average of k 2 cx t on a long, diffusive time scale, 
and thus Q oc i 2 . 

The flux and amplitude are obtained by assuming that 
the statistics of the turbulent flow v x are spatially ho- 
mogeneous and temporally stationary with the following 
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correlation function: 

(^(ki, ti)« x (ka, i 2 )> = (2^) 2 5(k! + k 2 )0(k 2 , t 2 - t t ) , 

(6) 

where <fi is the correlation function of v x in Fourier space. 
Further, the random turbulent flow is taken to have char- 
acteristic frequency u> and correlation time r c = 1/7. 
Specifically, tp is taken to have Lorentzian frequency spec- 
trum centered around to with width 7 as 0(k2, ti —t\) = 
V>(k2)/(^'/7r)e-^'( t2 - tl ) 7 /[(w'-w) 2 +7 2 ] with u > 7, 
where c.c. denotes complex conjugate. It can be read- 
ily shown that the velocity amplitude is related to the 
power spectrum ip as (v x ) = J d 2 kip(k)/ (2ir) 2 . Since 
k x ^ modes are generated by the shear [see, Eq. ©], ip 
is assumed to be dominated by modes with k x <C k y for 
simplicity. Then, the flux and amplitude of fluctuation 
of scalar fields, averaged over one oscillation of the shear 
flow (2tt/lu z ) as well as over the statistics of the turbu- 
lence follow from Eqs. Q and I© (see Eqs. (9)-(10) in 
@) and can be expressed in the following dimensionless 
form: 

( n ' V *)t=-J^ J d2fc ^(k) J J dTdTlA{P,W,T,Tl) 

(7) 

and 



Here. 

A(f3,TD, t, n) = cos{/?a(cos(r) — cos(ri)) — uj(t — n)}, 
and 



B(a,T,n) = (T-Ti){l + a 2 (I + icos2-ri) 

a 2 . 

+ — (sin 2r + 3 sin2ri — 8 cosri sinr)}. 

The dimensionless variables in Eqs. |(7]) and ([8]) are de- 
fined using wj 1 as a unit of time (i.e., r = w z i) as follows: 

- a = (x rms shearing rate 

- 7 = -2- = -^!— oc decorrclation rate of turbulence 

- j3 — kx oc scale separation between the mean and fluc- 

tuations 

- uj — — oc frequency of turbulence 



- to — -Tfifri oc molecular diffusion time scale 

Note that k = k y in is the typical wavenumber of 
turbulence in the y direction while k x (t) — k x (t\) + 
k y L dt'fl m sin(w z i') evolves in time due to shearing [see 
Eq. ©]. 

It is important to note that the flux can be expressed 
by using turbulent diffusivity Dt as {n'v x )t = —DtOxNq 
in homogeneous and stationary turbulence. In the case of 
inhomogeneous turbulence (e.g., due to the background 
density stratification or gradient in turbulence intensity) , 
non-diffusive flux can appear in the flux [for example see 
Ref [HI]. Further, Dt is likely to be positive, especially 
for a short-correlated turbulence or for weak turbulence 
quenched by strong shear, where a quasi-linear analysis 
becomes exact. Thus, in our case, the turbulent diffusiv- 
ity increases the dissipation rate of large-scale component 
Nq from the molecular value D to a larger value D + Dt- 
In the absence of shear flow, or equivalently, in the limit 
a — » 0, turbulent transport due to background turbu- 
lence is fast with the typical value of Dt ~ ul, where u 
and I are the characteristic velocity and length scale of 
turbulence. In the following sections, we will investigate 
how the turbulent transport and turbulence amplitude 
are affected by an oscillatory shear flow via numerical in- 
tegration of Eqs. and © by varying the values of a, 
7, (3, ZJ and td- In all cases, the molecular dissipation 
time is assumed to be large with tjj ^> 1. 



A. Short correlation time/period 

In the limit of short correlation time r c of turbulence 
(7 — > 00), double integrals of e~~'( T ~ Tl ' > in Eqs. (1) and 
(2) indicate that both the flux and the amplitude decay 
oc 7 _1 . Similarly, in the large frequency limit (uJ — > 00), 
both quantities scale as oc cJ~ 2 . In this regime, there is 
no possibility of resonance between shear flow and tur- 
bulence, with turbulence level and transport decreasing 
as either 7 or ZJ increases. This is because the time scales 
for the shear flow (t z = to' 1 and fi" 1 ) are far too large 
in comparison with the small correlation time (r c = 7 -1 ) 
or characteristic period (lu^ 1 ) of turbulence. Such a tur- 
bulent fluid is not able to recognize the subtle changes 
induced by the shear flow over its correlation time or 
period and is thus not affected by them. For instance, 
7>1 implies that 7 ^> uj z <==> r c <C r 2 , i.e., the tur- 
bulence decorrelates too rapidly to be influenced by the 
oscillation of the shear flow. 



T ^y- J d 2 kip(k) J J J dTdT 2 dTiA(/3,LJ,Ti,T 2 )UI. REGIMES OF DIFFERENT SCALINGS 

xe In this section, we identify the three regimes of differ- 

(8) ent scalings depending on the order of the relevant time 
scales. 



4 



B. Medium correlation time 

When turbulence has a longer correlation time with a 
finite value of 1 < 7 <C 00 in the limit of small dissipation 
t d > 1, the flux in Eq. ([7]) and amplitude in Eq. ([8]) 
take their maximum values when A{f3,W, r, n) = 1, i.e. 
when 



(3a [cos(t) — cos(ti)] — uj(t — n) = 0. (9) 

To obtain the condition for the occurrence of maxima, 
Eq. [5J we note that for 7 > 1, most of the contributions 
to the flux and amplitude in Eqs. ([7]) and ([5]) come from 
the time integral for |r— t\\ < 1. We can thus express cos 
in terms of sin in Eq. (O, expand sin as a Taylor series, 
and then factorize the first two leading order terms in the 
series for |t — ti| 1 to obtain the following condition: 

ZJ = ]3ac ^=4> uj/k-c^J(U%)t = 0. (10) 

Here, c = |/(t, ti) < 1 is a constant of less than unity, 
whose exact value depends on other parameter values. 
Since ui/k — uj/k y is proportional to the phase speed of 
turbulence, it is apparent that resonance occurs when the 
phase speed is approximately equal to the rms velocity of 
the shear flow. This is equivalent to the Doppler shifted 
frequency uid with the rms velocity of the shear flow being 
zero: 



Lu d ~ to - cky^J {U (x,t) 2 ) t = u - cf3Vt m w 0. 

(11) 

Note that a similar resonance condition in terms of rms 
velocity also holds in the case of random shear flow Q 
while it was overlooked in the analysis of an oscillatory 
shear flow @. In the case of a steady shear flow U^y, 
[H, the resonance condition pT|) becomes exact as aid — 
u> — kyllo = 0. 



C. Large correlation time 

When the turbulence correlation time is even longer 
with the value 7 <C 1, the flux and amplitude in Eqs. 
© and © can have considerable contributions from the 
time integrals for all values of \t—t\ | , leading to the possi- 
bility of the parametric type of resonance when ui = nui z 
for integer n. As indicated by Kim Q, the parametric 
resonance can be found by expanding e l acosT [ n terms 
of Bessel functions. We confirm numerically that Eqs. 
(|7|) and ([5]) do take maximum values when u> = nui for 
integer n in this case. 



1 


/ 


•< Enhanced Turbulent diffusion 


1 n 










3 0.8 




x v Background lurbulent dmusion- 


| 0.6 






D 






x 0.4 






0.2 






0.0 




Suppressed turbulent diffusion 



/^0 2 4 6 8 

/ a 

Analytically predicted Doppler resonance 

FIG. 1: Flux (dotted line) and amplitude (solid line) normal- 
ized to the a = case 
Figure 1 caption: Flux (dotted line) and amplitude (solid 
line) normalized to the a = case 



IV. NUMERICAL RESULTS 

In this section, we provide the quantitative predictions 
for turbulence level and turbulent transport in different 
regimes identified in §111, by numerical computations. 
The numerical accuracy was checked by benchmarking 
the code against an exactly integrable function which 
replicates the behaviour of our actual flux and ampli- 
tude. The numerical error was found to be within 0.0001 
and 0.001 for the flux and amplitude respectively. 

A. Medium correlation time 

We consider the medium correlation time in the limit 
of small molecular dissipation (large td ) . Specifically, we 
take the value 7 = 5 and td — 100. The scale separation 
is chosen to be f3 = 10. By using these parameter values 
and by varying a and 57 within the ranges of < a < 10 
and < ZU < 100, we perform numerical integrations of 
Eqs. ([7]) and ([5]) to obtain the scalings of turbulence am- 
plitude and transport with the shearing rate TO . These 
parameter values remain fixed unless stated otherwise. 

1. For a given Zo 

For a given turbulence frequency uJ, we normalize tur- 
bulence amplitude and transport by using values pro- 
vided by background turbulence in the absence of shear 
flow (i.e. for a = 0). Specifically, we choose the frequency 
ratio Zo = 10 and plot in Fig. 1 the flux and amplitude 
normalized to the a = case as a function of the rms 
shearing rate within the range of a £ [0, 10]. We note 
that as a consequence of our normalization, the value of 
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FIG. 2: Flux as a function of W — uj/ui z with a increasing 
from top to bottom: top line no shear, bottom line strong 
shear 



unity for the flux and amplitude at a = corresponds 
to that given by background turbulence in the absence 
of shear flow, with smaller or larger values describing 
the suppression or enhancement of turbulent diffusion, 
respectively. It is clear from Fig. 1 that as a increases 
from zero, the flux increases, indicating the enhancement 
relative to the case without a shear flow until it reaches 
its maximum, roughly at the Doppler resonance point, as 
noted in §111 B. As the shearing rate increases further be- 
yond this resonance point, the flux becomes quenched by 
strong shear. Similar behaviour, correlating well with the 
flux maximum, is also found in the amplitude although 
it is far less obvious. 



2. Flux for different ui 

In the absence of shear flow (a = 0), the flux due to 
background turbulence depends on the turbulence fre- 
quency. Thus, in order to examine the flux for differ- 
ent values of ZJ whilst varying a, we normalize the flux 
by its value in the case without a shear flow (a = 0) 
and for zero turbulence frequency ZJ = 0. With this 
normalization, we plot the flux as a function of the fre- 
quency ratio ZJ for different a S {0,0.2, ...,4} in Fig. 2. 
Fig. 3 shows the corresponding log plot. From Figs. 
2 and 3, one can see the two regimes of clearly differ- 
ent scalings. First, in the limit of large turbulence fre- 
quency, the flux always eventually decays oc jjj -1 - 96 ^- 06 . 
This is consistent with the analytically predicted behav- 
ior oc ZJ~ 2 in §111 A. On the other hand, the maximum 
value of the flux at the Doppler resonance point decays 
oc £j-°- 486 ±o.oi^ much glower with increasing ZJ. Fig. 4 
shows similar behavior as rms shearing rate a increases. 
Specifically, in the strong shear limit, flux eventually de- 
cays proportionally to a -i-0i5±o.oi ^ Q-i.oisio.oi^ m 

agreement with Kim 0] while the resonant flux decays 
proportionally to a -o.534±o.oi K Q-o.534io.oi_ There _ 



FIG. 3: Logarithmic plot of Fig [2] 
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FIG. 4: Flux as a function of a for T3 = ll>/lj z with 57 in- 
creasing from top (random turbulence) to bottom (wave like 
turbulence) 

fore, the scaling of the flux with Q m is weaker for max- 
imum resonance flux with roughly —1/2 power-law de- 
pendence. Furthermore, the absolute value of the flux at 
the resonance is larger than that due to background tur- 
bulence, manifesting the resonance enhanced transport. 
In other words, oscillatory shear flow can enhance the 
overall transport above the value given by background 
turbulence due to resonance while it quenches transport 
for sufficiently strong rms shear. The resonance enhanced 
transport becomes more apparent if the flux is normal- 
ized to the case a — case for all possible values of 
turbulence frequency uj in which case turbulent diffusion 
scales as oc a 1 ' 420 * ' 1 ^ Qi.420io.i_ That ig _ turbulent 
diffusion increases because of the oscillatory shear flow! 



3. Amplitude for different u> 

Fig. 5 shows turbulence amplitude as a function of 
ZJ £ [0,40] for different values of a <E {0,0.2, ...,4}, again 
normalized to the a — to — case. In the limit of 
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FIG. 5: Turbulence amplitude as a function of u> = w/o; z with 
a increasing from top (no shear) to bottom (strong shear) 



FIG. 7: Relationship between [a, ui] for resonant points with 
/3 increasing from top to bottom 
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FIG. 6: Turbulence amplitude as a function of a for different 
values of T3 with Zj increasing from (random turbulence) to 
bottom (wave like turbulence) 



large turbulence frequency, the amplitude decreases as 
oc (j- 1 - 901 ^ - 1 while the resonant amplitude decays less 
rapidly with u; as oc £j-°- 476±0 - 05 . On the other hand, Fig. 
6 is a log plot of the amplitude as a function of a, showing 
the two regimes with clear scalings with a. In the strong 
shear limit, amplitude oc a -i-402±o.05 K q-1.402±o.05 
while the resonant amplitude oc a -o.533±o.05^-o.533±o.o5 
Therefore, maximum amplitude has a slower decrease as 
shearing rate increases compared to asymptotic value of 
amplitude in the limit of strong shear. 



4- Cross-phase results 



malized flux (the so-called cross-phase) cos 9 as 

(n'v x ) 



cos V = 

By using the results obtained in the previous subsections, 
we obtain the scaling of cross-phase with uJ in the large 
frequency ZJ limit and at the Doppler resonance points 
as: 



cost 
cost' 



OC 0J 



-1.009±0.035 



— -0.248±0.035 
OC LO , 



respectively. Similarly, in the strong shear limit and at 
the resonance points, the cross-phases scale with a (oc 



cost 
cos( 



oc Q 
oc 



-0.315±0.035 



q-0.267±0.11 



(12) 



respectively. The integration accuracy is always less than 
0.0001 and 0.001 for the flux and amplitude accordingly, 
consequently errors in the scalings were found and related 
to the cross-phase. 

In comparison with a steady shear flow case where 
cos6> oc ft -1 / 6 HIH, cos6» in Eq. (JUJ) have stronger de- 
pendence on £l m , suggesting that oscillatory shear flow 
is more efficient than a steady shear flow at reducing 
normalized transport by affecting the phase relation be- 
tween n' and v x . In particular, cost9 oc J7^°- 315±0 035 in 
the strong shear limit, with its value decreasing rapidly 
compared to the steady shear flow case [4j. Therefore, 
oscillatory shear flows are more efficient in reducing nor- 
malized flux by quenching flux more than amplitude. 



The flux (n'v x ) involves not only the phase relation 
between n' and v x but also turbulence amplitude. We 
quantify the flux due to phase shift between n' and v x 
only, independent of turbulence level, by defining the nor- 



5. Behaviour of the Doppler resonance 

Figure 6 shows us the relation between the rms shear- 
ing rate of the oscillatory shear flow and turbulence fre- 
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FIG. 8: Flux as a function of u for ui z — 0.5 (dotted line) 
and lo z = 0.125 (solid line) when 7 = 0.01, showing multiple 
harmonic resonance in addition to Doppler resonance 



FIG. 9: Flux as a function of lo for lo z = 0.5 (dotted line) 
and 10 z — 0.125 (solid line) when 7 = 10~ 4 , showing more 
pronounced resonant peaks than Fig. 8 



quency at the resonance for different values of the scale 
separation. In each case the relationship appears to be 
linear, with its proportionality decreasing as scale sep- 
aration increases. Specifically, the relationship between 
these quantities at the peak resonance is found to be 

_ 7r-l.OOliO.Ol . „ _ ,_ -5-0.944±0.05 „ 

lo oc p a lor tnc Hux and lo oc p a lor 

the amplitude. This is consistent with analytical esti- 
mate in §111 within the numerical error caused by finite 
computing time. 



B. Long correlation time limit 

For sufficiently long correlation time of turbulence, 
parametric type resonance between oscillatory shear flow 
and turbulence can occur where lo = mo z for integer n, 
as discussed previously in §111 C (see also [9(). In order 
to explore this resonance, we choose sufficiently small 7 
with a specific value of 7 = 0.01 whilst taking the dissi- 
pation time to be td = 10 3 . Furthermore, we change our 
non dimensional variables so that lo and lo z are the only 
free parameters and alter the time average from being 
over [0,2tt] to [0,2tt/w z ]. 

Figure 8 is the plot of the flux as a function of lo for 
lo z = 0.5 (dotted line) and lo z = 0.125 (solid line) for the 
fixed values 7 = 0.01 and td = 10 3 , and shows that the 
resonance peaks are equally spaced and increase in inten- 
sity until it reaches the transition phase into the Doppler 
resonance. Interestingly, the flux at Doppler resonance 
exceeds that due to parametric resonances. Beyond the 
Doppler resonance point, the flux decays rapidly to zero 
as shearing rate increases. Figure 9 shows how reduc- 
ing the decorrelation rate (to 7 = 10~ 4 ) and molecular 
dissipation (to to = 10 4 ) acts to amplify the height of 
these peaks as further harmonics appear. The height of 
resonance peaks would blow up in the limit of infinite 
memory time of turbulence (i.e. no stochasticity 7^0) 
and no molecular dissipation (i.e. td — > 00), leading to 



parametric instability [16[ . Parametric resonance can be 
physically understood since the equilibrium provided by 
a large-scale shear flow with frequency lo z requires fluc- 
tuations to be invariant under the time translation by 
1/lo z , thereby supporting the excitation of modes of fre- 
quencies of lo = ruo z . It is worth noting that a classical 
example of parametric resonance is a vertically oscillat- 
ing pendulum [2(| • The pendulum sways side to side and 
the frequency of these oscillations are integer multiples 
of the driving frequency, depending on the magnitude of 
the oscillations. For instance, the case when the driving 
frequency and the pendulum oscillation frequency match 
is simply explained by the driving frequency trying to 
force the pendulum into step. 

V. CONCLUSION 

We have performed a detailed study of the effect of os- 
cillatory shear flow on turbulent transport and amplitude 
in passive scalar field model. Specifically, we have identi- 
fied the three distinct scaling regimes with rms shearing 
rate, depending on the value of turbulence decorrelation 
time. 

• Short correlation time limit (7 — > 00): there is no 
resonance due to too rapid change in turbulence 
characteristics, with the flux decaying as 7 -1 . 

• Medium correlation time (1 < 7 <C 00): turbulent 
transport is enhanced around the rms Doppler reso- 
nance where the wave phase speed matches the rms 
velocity of oscillatory shear. Increasing either the 
shear or the characteristic frequency from Doppler 
resonance point acts to quench turbulent transport. 
Our numerical results with 7 = 5 show that if 
the flux is normalized to the case of a = and 
lo = 0, the flux, thereby turbulent diffusion, scales 
as, oc f2~ 05343 at resonance points. The depen- 
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dence on f2 m at resonance point is thus weaker than 
that in the strong shear limit (oc Sl~ 1015 ). Further- 
more, the value of flux at resonance point is larger 
than that in the absence of the shear flow, mani- 
festing the enhancement of turbulent transport due 
to the oscillating shear flow. If the flux is normal- 
ized to the case a = 0, the flux scales as oc ^„ 831 
for all possible values of turbulence frequency w, 
highlighting the increase in turbulent diffusion due 
to oscillatory shear flow. 

• Long correlation time limit (0 < 7 < 1): there 
is parametric (harmonic) resonance of the form 
lu = nuj z where n 6 N in addition to Doppler reso- 
nance. At these resonance points, the numerically 
computed flux and amplitude show significant en- 
hancements, with their maximum values provided 
by Doppler resonance peaks. 

These results are summarised in Table I. Our results 
suggest that oscillatory shear flows can enhance the value 
of turbulent diffusion, speeding up turbulent mixing, de- 
pending on the characteristics of shear flow and turbu- 
lence. Enhancement of turbulent mixing can then have 
either welcoming or unwelcoming consequences in turbu- 
lent mixing. A typical example of the latter can be found 
in laboratory plasmas where the confinement is a critical 
issue. These results thus suggest that for the understand- 
ing and predictive modelling of turbulent transport in 
plasmas, it is necessary to determine frequency and power 
spectra of shear flows and turbulence. Our results have 



implications for turbulent mixing in many other fields 
such as geophysics, oceanography, atmospheric physics, 
solar physics, and magnetohydrodynamics where shear 
flows and turbulence are main players in transport. In 
particular, similar results obtained in the passive scalar 
fields model are expected to be valid for the transport 
of magnetic fields in the 2D magnetohydrodynamic tur- 
bulence as long as the backreaction of magnetic fields is 
negligible (i.e. in the kinematic regime). An interesting 
question is then what happens to the transport / diffusion 
of large-scale magnetic fields when the backreaction is 
sufficiently strong to modify the characteristics of tur- 
bulence. In particular, it would be interesting to study 
whether oscillatory shear flows can weaken the severe 
quenching in the amplification of magnetic fields (the so- 
called dynamos) in 3D and their diffusion rate in 2D due 
to the magnetic backreactions. Furthermore, the back- 
reactions of turbulence on shear flows will also have an 
importance consequence on the evolution of shear flows, 
for instance, leading to complex temporal and spatial 
dynamics [131 ], thereby dynamically determining the fre- 
quency and power spectra of shear flows, which have been 
assumed to be given in this paper and the previous works 
[H H 0, Is! ■ These problems will be addressed in future 
publications. 
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7 — > oo No resonance 


(n' 2 ) oc uj~' a , cc 7 _i 
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Doppler 
resonance 
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TABLE I: Summary of results with < 7 < 1, 1 < 7 « oo 
and 7 — > oo describing the long, medium and short correlation 
times respectively. 



